Mechanism of hyperscaling violation in the 2D 8-state Potts model with long-range 
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The first-order phase transition of the two-dimensional eight-state Potts modei is shown to be 
rounded when long-range correlated disorder is coupled to energy density. Critical exponents are 
estimated by means of large-scale Monte Carlo simulations. In contrast to uncorrelated disorder, a 
violation of the hyperscaling relation j/v — d — 2x„ is observed. This violation is caused by large 
disorder fluctuations, like in the 3D random field Ising model. In the thermal sector too, evidences 
are given for such violation in the two hyperscaling relations a/v = d — 2x e and Ijv — d — x t . The 
scaling dimension of energy is conjectured to be x € = a/2, where a is the exponent of the algebraic 
decay of disorder correlations. 
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Quenched disorder when coupled to the energy den- 
sity, say by dilution or random couplings, is known to 
soften first-order phase transitions. As argued by Imry 
and Wortis local fluctuations of impurity concentra- 
tion can destabilize the ordered phases in coexistence at 
the transition temperature if the surface tension is suf- 
ficiently small. In 2D, it was rigorously proved that an 
infinitesimal amount of disorder is sufficient to make any 
first-order transition continuous [2j,|3j]. The complete van- 
ishing of the latent heat was first observed numerically 
in the case of the 2D 8-state Potts model The criti- 
cal behavior of the disorder-induced second-order phase 
transition is governed by a new random fixed point [Bj]. 
The universality class was later shown to depend on the 
number of states q 0. In 3D, a finite disorder is required 
to round completely the first-order phase transition. The 
phase diagram exhibits a tricritical point separating a 
first-order regime from the disorder-induced continuous 
one, as first observed in the bond-diluted 4-state Potts 
model • A rounding of the first-order phase transition 
of the 2D Potts model was also reported for anisotropic 
aperiodic sequences of couplings [8J , and for layered ran- 
dom couplings 0, In both cases, the couplings are 
infinitely correlated in one direction. In the random case, 
the critical behavior was shown to be governed by a q- 
independent infinite-randomness fixed point. The critical 
exponents are therefore those of the layered random Ising 
model, the celebrated McCoy- Wu model llj, ll2| . Inter- 



estingly, the same critical behavior is observed for the 
Potts model with homogeneous uncorrelated disorder in 
the limit q +oo 13;] ■ Hyperscaling holds for all these 
models. 

In the following, we consider random bond cou- 
plings Jij > with algebraically decaying correla- 
tions (J(0) — J)(J(r) — J) ~ r~ a . A simple general- 
ization of the Imry- Wortis criterion shows that the low- 
temperature phase is destabilized when the coupling fluc- 



tuations, inside a domain of characteristic length £, 
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£ d ~ a/2 , (a < d) (1) 



increase faster with t than the interface free energy at d ~ Y . 
In the two-dimensional case, we expect the first-order 
phase transition to be softened for a < 2. 

We consider the 2D g-state Potts model with Hamil- 
tonian 



ihi) 



(2) 



where <jj £ {0,1, . . . ,q — 1} and the sum extends over 
pairs of nearest neighbors of the square lattice. We re- 
strict ourselves to the case q — 8 for which the correlation 
length of the pure model is £ ~ 24 at the transition tem- 
perature. We considered a binary distribution of coupling 
constants Jij 6 {Ji, J%] with 

(e Jl -l)(e*-l)=«. (3) 

In the case of uncorrelated disorder, Eq. ^ is the self- 
duality condition that gives the location of the critical 
line. The ratio r = J%j J\ is used as a measure of the 
strength of disorder. Here, we presents results for the 
case r = 8. To generate correlated coupling configura- 
tions {Jij}, we simulate another spin model, namely the 
Ashkin- Teller model (o"i,Ti = ±1) 
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at different points of its critical line e~ 2K T = sinh2J AT . 
Two symmetries of the Hamiltonian are spontaneously 
broken at low temperatures: the global reversal of the 
spins <Ji and the reversal of both at and t;. There- 
fore, two order parameters can be defined, magnetiza- 
tion Oi and polarization (JiTi, leading to two inde- 
pendent scaling dimensions /3 AT = 2 4-i6y anc ^ Par = 



2 



wherein we use the parametrization cos W- 



ponent is v 



1 (j £ [0;4/3]). The correlation length ex- 

AT = z-'iy • Spin configurations of this model 
are generated by Monte Carlo simulation using a cluster 
algorithm 14|. For each of them, a coupling configura- 
tion of the Potts model is constructed as 



Jl + Jl , Jl ~ J2 



-CTiTi, 



(5) 



where the site j is either at the right or below the site 
i. By construction, disorder fluctuations are self-similar 
and the coupling constants display algebraic correlations 
(Jij — J)(Jki — J) ~ \n — rk\~ a at large distances with 
a = 2/3^F ' /v AT — f/(4 — 2y). We have considered six 
points on the critical line, y € {0, 0.25, 0.50, 0.75, 1, 1.25}, 
leading to six correlated disorder distributions with a ~ 
0.25, 0.286, 0.333, 0.4, 0.5 and 0.667. Note that for y = 0, 
the Ashkin- Teller model is equivalent to the 4-state Potts 
model. The critical behavior is therefore affected by log- 
arithmic corrections. Finally note that this construction 
ensures that the constraint ([3]) implies the self-duality of 
our random Potts model. 

The Potts model is then simulated using the Swendsen- 
Wang algorithm 15j. Lattice sizes between L = 16 and 
256 are considered. For each disorder configuration, 1000 
MCS are performed to thermalize the system and 20, 000 
MCS for data accumulation (auto-correlation time is 
r ~ 2 for L = 256). Thermodynamic quantities are av- 
eraged over a number of disorder configurations propor- 
tional to 1/L 2 . For the largest lattice size L = 256, 2560 
disorder configurations are generated while for L = 64 
for instance, this number is raised up to 40960. Stabil- 
ity of disorder averages is checked. In the following, we 
will denote (X) the average of an observable over ther- 
mal fluctuations and (X) the average of the latter over 
disorder. 

On the critical line, the typical spin configurations of 
the Ashkin- Teller model display a large cluster of polar- 
ization err = +1 or —1. As a consequence, our random 
coupling configurations also exhibit large clusters of ei- 
ther strong or weak bonds. The probability distribution 
of the total energy of the Potts model shows two peaks 
corresponding to these two kinds of bond configurations. 
This distribution is highly correlated to the probability 
distribution of polarization of the Ashkin- Teller model. 
Since the latter undergoes a second-order phase transi- 
tion, the two peaks come closer as the lattice size is in- 
creased. Note that in our model, macroscopic region of 
strong couplings are not rare: they have a probability 
1/2. Moreover, they have a fractal dimension 1 < dj < 2 
determined by the Ashkin- Teller model and, in the ther- 
modynamic limit, only one such macroscopic region is 
expected to be present in the system. For this reason, 
the transition is not smeared but rounded (lif . 

We estimate critical exponents by Finite-Size Scaling. 



The exponent (3/v can be extracted from magnetization 
(m) and its moments (m n ) with n — 2,3,4. We observe 
nice power laws without any significant correction to scal- 
ing. Our estimates of x a = fijv evolve with a and range 
from 0.061(5) (y = 0) to 0.108(4) (y = 1.25), to be com- 
pared with 0.150(2) for uncorrelated disorder (Tab. IJ). 
We then consider the average magnetic susceptibility, nu- 
merically computed via the fluctuation-dissipation theo- 
rem x = L d (m?) — (to) 2 . The data display large cor- 
rections to scaling (see Fig. [1]). A cross-over is observed 
around L = 48, not far from the correlation length £ = 24 
of the pure 8-state Potts model. In the region L > 96, 
power-law interpolations give stable estimates for j/v 
going from 1.70(7) (y = 0.00) to 1.62(4) (y = 1.25), 
to be compared with 1.69(4) for uncorrelated disorder. 
The hyperscaling relation 7/1/ = d — 2j3/v is therefore 
not satisfied. We shall identify disorder fluctuations as 
the origin of this hyperscaling violation, like in the 3D 
Random-Field Ising Model (RFIM). Consider the follow- 
ing decomposition: 



X = L d [(m 2 ) - (m) ] - L d [(m) 2 - (m)' 



(6) 



The first term, when computed separately, displays a 
power law behavior with an exponent (7/^)* incompat- 
ible with 7/1/ but in agreement with the hyperscaling 
relation (Tab. H]). The second term of ©, the so-called 
disconnected susceptibility, involves the ratio 



Rr, 



(m) 2 — (to) 
(to) 2 



(7) 



which is expected to behave as R m ~ R m (oo) + AL^^, 

if magnetization is not self-averaging [l^. We indeed 

observe that R m goes to a non-vanishing constant in the 

limit L — > +00 (Fig. [5]). The second term of Eq. © 
2 

therefore behaves as L d {m) ~ L d - 2 P/ V j i. e . with an 
exponent satisfying the hyperscaling relation. Since the 
data indicate that it is also the case for the first term 
of Eq. (|6]), one can imagine that, if their amplitudes 
are equal, the dominant terms will cancel. To test this 
hypothesis, we compute the ratio 



(to 2 ) — (to) 



(to) 2 — (to) 



(8) 



We observe a plateau at 1.00(4) (Fig. [5]). We therefore 
conclude that the hyperscaling violation is the result of 
an exact cancellation of the dominant contributions of the 
two terms of Eq. ©. The scaling of the average suscep- 
tibility is therefore determined by the first non- vanishing 
scaling correction % ~ L d ~ 2 ^l v ^ u) of any of the two terms 
of Eq. ((6]) and the hyperscaling violation exponent 9 is 
the exponent to of this correction. Note that this is also 
the mechanism of hyperscaling violation invoked in the 
context of the RFIM [TJ| where our exponent (7/f)* is 
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FIG. 1. Average magnetic susceptibility \ (bottom) and 

L d [(m 2 } — (m) ] (top) versus lattice size L for different values 
of y and uncorrelated disorder (Uncorr.). The dashed line 
indicates that the fit is performed over lattice sizes L > 96 
only. 



as 



0.6 

I 

|A 

'.vo.4 



0.2 




- 1.6 






v=0.00 


□ 


y=0.25 





V=0.50 




y=0.75 


< 


y=1.00 


V 


y=1.25 




Uncorr 



- 1.4 



V 

I 

"'A 
V 

-I 
V 

I 

A 



0.02 0.04 
1/L 



** * * 

i i i 

0.02 0.04 
1/L 



- 1 



FIG. 2. Ratios defined by Eq. (J7J) on the left and 
right versus the inverse of the lattice size. 
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denoted 4 — fj 19]. In the case of uncorrelated disorder, 
the ratio goes to a value significantly different from 
1 in the large size limit (Fig. [2]). The dominant contri- 
bution of the two terms of ([6]) do not cancel in this case 
and therefore hyperscaling is not violated. 

The divergence of specific heat is completely washed 
out by the introduction of disorder, which means that 
the specific heat exponent a/v is zero or negative (Fig[3]). 
C = L d < e 2 > — < e > 2 can be decomposed in the same 
way as %: 



C = L d [(e 2 ) - (e) ] -L d [(e)i-(ey 



(9) 



We observe a nice power-law behavior of the first term 
with an exponent in good agreement with (a/v)* — 
(j/v) AT = d — a, which means that the fluctuations of 
energy are dominated by the fluctuations of the couplings 
and therefore of the polarization density in the original 
Ashkin- Teller model. The second term involves the ra- 
tio R e , constructed in the same way as R m Q. Our 
numerical data show that energy is not self-averaging 



FIG. 3. Average specific heat C (bottom) and L d [(e 2 ) - (e) ] 
(top) versus lattice size L for different values of y and uncor- 
related disorder (Uncorr.). For clarity, error bars of C in the 
uncorrelated case have been drawn as dashed line when they 
overlap with other points. 
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FIG. 4. Ratios defined in the same way as Eq. (0 on the left 
and on the right but for energy instead of magnetization. 
On the right, the error bars for uncorrelated disorder are large 
(they overlap other bars for large lattice sizes) and have not 
been represented for clarity. 



(i? e (oo) ^ 0) and the ratio 



<e z >-<e>- 



exhibits a plateau 



<e>^-<e>" 

at the value 1.00(5) (Fig. HJ. This implies the cancella- 
tion of the dominant contribution of the two terms of C so 
that a violation of the hyperscaling relation ajv = d—2x e 
is expected. Even though we cannot measure x e from 
the scaling behavior of energy, we infer that it can be ex- 
tracted from the hyperscaling relation (a/v)* = d — 2x e 
which implies x t = a/2. In the case of uncorrelated disor- 
der, R e is compatible with zero which means that energy 
is self-averaging and therefore the two dominant contri- 
butions of Eq. ([9]) do not cancel. As observed, hyper- 
scaling is not violated in this case. 

In pure systems, a good estimator for the deter- 
mination of the correlation exponent v is — d h ^ ? "^ — 



- d (me)-(m)(e) 



This is generalized to random systems 
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FIG. 5. Quantity 



d In (m) 
dp 



(bottom) and L 



d (me) — (m) (e) 
7m7 



(top) versus lattice size L for different values of y and un- 
corrected disorder (Uncorr.). 
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FIG. 6. Ratios defined by Eq. (JT2} on the left and (JT3} on the 
right. On the right, the error bars for uncorrelated disorder 
are large and have not been represented for clarity. 
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and is expected to scale as d — x e . A power-law interpo- 
lation of our data gives exponents 1/v close to zero but 
with large error bars. Consider again the decomposition 



d In (m) 

dp 



= L° 



(me) - (to) (e)j , [(m)(e) - (to) (e) 



-L 



(to) 



(11) 

The first term displays a power-law behavior with expo- 
nents 1/v* close to, though slightly above, d — x € , which 
implies x e ~ a/2 (Tab. fl}. This estimate is consistent 
with the one obtained from the specific heat. The sec- 
ond term of (fTTj) involves the ratio 



Rr, 



(m)(e) — (to) (e) 

MM ' 

,(L) ~ R me (oo) + aL^' with 



(12) 



which behaves as R n 
0.3. The constant R me (oo) is clearly finite, except maybe 
for y — 1.25. Like in the magnetic case, the two terms of 
Eq. (ITT1) have the same dominant scaling behavior. The 
ratio 



(me) — (to) (e) 



(to) (e) - (to) (e) 



(13) 



displays a plateau at 1.00(4) (Fig. [5]). Consequently, the 
dominant contribution of the two terms of (jlll) is the 
same and they cancel. Hence, the hyperscaling relation 
\/i> = d — x e is expected to be violated. In the case of 
uncorrelated disorder, R me is compatible with zero so we 
do not expect any cancellation of the two dominant con- 
tributions of Eq. (|11[) and, consequently, no hyperscaling 
violation. 



y 


0.25 0.5 0.75 1 1.25 


a 


0.25 0.286 0.333 0.4 0.5 0.667 


d-2p/v 

(jM* 


0.061(5) 0.060(5) 0.067(5) 0.075(5) 0.091(5) 0.108(4) 
1.88(1) 1.88(1) 1.87(1) 1.85(1) 1.82(1) 1.784(8) 
1.70(7) 1.69(8) 1.69(7) 1.67(8) 1.66(6) 1.62(5) 
1.91(2) 1.90(3) 1.89(3) 1.87(3) 1.83(3) 1.79(3) 


(a/vf 
d — a 


1.75(1) 1.73(2) 1.68(2) 1.61(2) 1.51(2) 1.34(2) 
1.75 1.714 1.667 1.600 1.500 1.333 


l/v* 
d-a/2 


1.90(2) 1.89(2) 1.86(2) 1.83(2) 1.78(2) 1.69(2) 
1.875 1.857 1.835 1.8 1.75 1.667 



TABLE I. Critical exponents measured by Monte Carlo sim- 
ulations, or computed from them, and conjectured values. 

Finally, note that, in contrast to the RFIM, at least 
two hyperscaling violation exponents are needed to ex- 
plain the numerical data presented in this letter. In the 
magnetic sector, the hyperscaling violation exponent is 
estimated to be 8 = (7/^)* — {l/ v ) — 0.2 while in the 
energy sector, 9 — (a/v)* — {a/v) >d — a. 
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